Highly Damped Quasinormal Modes of Kerr Black Holes: 
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We compute for the first time very highly damped quasinormal modes of the (rotating) Kerr 
black hole. Our numerical technique is based on a decoupling of the radial and angular equations, 
performed using a large- frequency expansion for the angular separation constant s Ai m . This allows 
us to go much further in overtone number than ever before. We find that the real part of the 
quasinormal frequencies approaches a non-zero constant value which does not depend on the spin s 
of the perturbing field and on the angular index I: ujr = mzu(a). We numerically compute -cu(a). 
Leading-order corrections to the asymptotic frequency are likely to be ~ l/u)i. The imaginary part 
grows without bound, the spacing between consecutive modes being a monotonic function of a. 
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I. INTRODUCTION 

Black holes (BHs), as many other objects, have char- 
acteristic vibration modes, called quasinormal modes 
(QNMs). The associated complex quasinormal frequen- 
cies (QN frequencies) depend only on the BH fundamen- 
tal parameters: mass, charge and angular momentum. 
QNMs are excited by BH perturbations (as induced, for 
example, by infalling matter). The early evolution of a 
generic perturbation can be described as a superposition 
of QNMs, and the characteristics of gravitational radi- 
ation emitted by BHs are intimately connected to their 
QNM spectrum. One may in fact infer the BH parame- 
ters by observing the gravitational wave signal impinging 
on the detectors |l| : this makes QNMs highly relevant in 
the newly born gravitational wave astronomy 0, • 

Besides this "classical" context, QNMs may find a very 
important place in the realm of a quantum theory of grav- 
ity. General semi-classical arguments suggest Q that on 
quantizing the BH area one gets an evenly spaced spec- 
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trum of the form 

An = 4 log (fc)/£n ; n = 0,1,... (1) 

where lp is the Planck length, and k is an integer to be 
determined. Hod proposed to fix the value of k, and 
therefore the area spectrum, by promoting QN frequen- 
cies with a very large imaginary part to a special position: 
they should bridge the gap between classical and quan- 
tum transitions. Hod obtained, for the Schwarzschild 
BH, k = 3. Following his proposal, further enhanced by 
the prospect of using similar reasoning in Loop Quantum 
Gravity to fix the Barbero-Immirzi parameter |fj , the in- 
terest in highly damped BH QNMs has grown consider- 
ably 0- There is now reason to believe that the connec- 
tion between QN frequencies and the BH area quantum 
is not as straightforward as initially suggested. How- 
ever a relation between classical and quantum BH prop- 
erties does seem to exist, even in non-asymptotically flat 
spacetimes @ • A prerequisite to study this connection is 
to compute QN frequencies having very large imaginary 
part. So far this problem has been solved only for a few 
special geometries: Schwarzschild BHs 0, 
Reissner-Nordstrom (RN)BHs [T1 111 liafthe Banados- 
Teitelboim-Zanelli BH [jjj, and the four-dimensional 
Schwarzschild-anti-de Sitter BH [l5| . 

We must try to include the Kerr geometry in this short 
catalogue, due to its great importance and simplicity. 
This is a problem of great relevance for the scientific 
community, and quite a lot of effort is being invested 
here. This effort is in direct proportion to the difficulty 
of the problem. All previous attempts 0, ITfl IltI ITsj 
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at probing the asymptotic QNMs of Kerr BHs have been 
unsuccessful, or at least unsatisfactory. There have been 
several contradictory "analytical" results, which were ei- 
ther based on incorrect assumptions, or could not probe 
the highly damped regime |16(. The few numerical re- 
sults 0> EH are not decisive either, although they 
definitely show some trend. A numerical investigation is 
necessary both as a benchmark and as a guide to ana- 
lytical approaches. Here we carry out such a numerical 
study. We improve on previous results by going further 
in overtone number than ever before, in order to really 
probe the asymptotic re gim e. The starting point for our 
analysis is, as previously [lall8l| . Leaver's continued frac- 
tion technique as improved by Nollert , with a few 
appropriate modifications |l3j |. However, we now decou- 
ple the angular and radial equations. We first determine 
the asymptotic expansion for the angular separation con- 
stant, and then replace this asymptotic expansion in the 
radial equations. This trick spares us the need to solve 
simultaneously the two equations, which was the major 
drawback of previous numerical works. A leading-order 
asymptotic expansion of the separation constant is typi- 
cally accurate for \aui\ > 5, where a is the dimensionless 
Kerr rotation parameter |2(| ■ In this study we go well be- 
yond this regime (we can usually compute modes up to 
| aw | > 50, an order of magnitude larger). So we have 
great confidence that our results really yield "asymp- 
totic" QN frequencies. 

We find that our previous results for negative m 
and moderately damped QN frequencies were quite close 
to the true asymptotic behaviour (especially for large val- 
ues of a) , while convergence to the asymptotic value was 
not yet achieved for positive m. Our improved calcu- 
lations have been carried out with two independent nu- 
merical codes. Our main results are that: (i) The real 
part of the QN frequencies lur approaches a non-zero 
constant value. This value does not depend on the spin 
s of the perturbing field and on the angular index I. It 
only depends on the rotation parameter a, and is pro- 
portional to m: lur ; = mm (a). We determine zu(a) nu- 
merically. A fit of our numerical data by power series in 
l/u>i and wl/ujj suggests that leading-order corrections 
to the asymptotic frequency should be of order (ii) 
The imaginary part u>j grows without bound, the spac- 
ing between modes 5u>i being a monotonically increasing 
function of a. 



II. BASIC EQUATIONS 

In the Kerr geometry, the condition that a given fre- 
quency be a QN frequency can be converted into a state- 
ment about continued fractions, which is rather easy to 
implement numerically. Such a procedure has been ex- 
plained thoroughly by Leaver [l9|. so here we shall only 
recall the basic idea. The perturbation problem reduces 
to a pair of coupled differential equations: one for the 
angular part of the perturbations, and the other for the 



radial part. In Boyer-Lindquist coordinates, defining 
u = cos 9, the angular equation reads 



[(l-u 2 )S lm , u ] u + 
(au>u) 2 — 2aujsu + s + s Ai„ 



(m + suY 
l-u 2 



(2) 



5*;™ — 0, 



and the radial one is 



AR (s + l)(2r - l)R lm , r + V{r)R lm = 0, (3) 



where A = r 2 — r + a 2 and 



V(r) 



(r 2 + a 2 ) 2 uj 2 + is(am(2r - 1) - uj(r 2 - a 2 )) + 



+a 2 m 2 — lamijjr 



A 1 + [Zisivr - a 2 oj 2 - s A lm ] . (4) 



The parameter s = 0,-1,-2 for scalar, electromagnetic 
and gravitational perturbations respectively, and s Ai m 
is an angular separation constant. In the Schwarzschild 
limit the angular separation constant can be determined 
analytically: s Ai m — 1(1 + 1) — s(s + 1). 

Boundary conditions for the two equations can be cast 
as a couple of three-term continued fraction relations |19|| . 
Finding QN frequencies is a two-step procedure: for as- 
signed values of s,l,m,a and u, first find the angular 
separation constant s Ai m (uj) looking for zeros of the an- 
gular continued fraction; then replace the corresponding 
eigenvalue into the radial continued fraction, and look 
for its zeros as a function of u>. This has been the strat- 
egy adopted in earlier works 0, 0, 0, , where the 
first ~ 50 modes were computed. These numerical inves- 
tigations showed a rich (and perhaps confusing) behav- 
ior. For negative to and large enough a the first ~ 50 
modes display some kind of convergence, and are consis- 
tent with our new results. Among modes with positive 
to, only those having \s\ = I = m = 2 seemed to con- 
verge. This convergence was deceiving: positive-TO re- 
sults in |l§| were not yet in the asymptotic regime. The 
"true" asymptotic behavior turns out to be much simpler 
than the intermediate-damping regime explored in |l8| . 

Since the major numerical difficulty lies in the coupling 
of the two continued fractions, here we adopt a "trick" 
to decouple them. We first carry out a careful study of 
the angular equation, determining the asymptotic form 
of the separation constant s Ai m for frequencies oj with 
large imaginary part. Then we substitute this expansion 
in Eqs. By this trick we reduce the problem to 

the numerical solution of a single three-term recursion re- 
lation: it is then possible to probe the asymptotic regime 
for highly damped modes. In the following section we 
shall briefly discuss our main analytical and numerical 
results for the asymptotic expansion of s Ai m . 
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FIG. 1: Real part of high-order QN frequencies for scalar 
(s = 0), electromagnetic (s = —1) and gravitational (s = —2) 
perturbations of a Kerr BH with a = 0.1 (I = m = 2) as a 
function of their imaginary part. QN frequencies of different 
spins converge to the same value. For any kind of perturbation 
we are already deep in the region of validity of the asymptotic 
expansion J7J, since \au\ ~ 60. 



FIG. 2: Scalar QN frequencies of a Kerr BH with a — 0.1, 
I = 2 and different values of m. The asymptotic value is 
proportional to m. Modes with m — oscillate around zero. 
The amplitude of these oscillations at fixed a decreases very 
fast, probably tending to zero as uii — > oo: this is consistent 
with the behavior shown in Fig. 6 of fTs|l . 



III. ASYMPTOTIC EXPANSION OF THE 
ANGULAR SEPARATION CONSTANT 

The analytical properties of the angular equation (J2J 
and of its ei genv alues have been studied by many au- 
thors [2ll2aMl2ll25|. Series expansions of s A hn for 
I aw I <C 1 are available, and they agree well with nu- 
merical results 25]. On the other hand, the asymptotic 
behavior for large frequencies has hardly been studied at 
all. An analytical power-series expansion for large (pure 
real and pure imaginary) values of aui can be found in 
Flammer's book |21| . but it is limited to the case s = 0. 
Flammer's results are in good agreement with the ex- 
haustive numerical work by Oguchi |26|. who computed 
angular eigenvalues for complex values of auo and s = 0. 
A review of numerical methods to compute eigenvalues 
and eigenfunctions for s = can be found in [27| . Quite 
surprisingly, there are no systematic numerical results 
for general spin s, and the few analytical predictions for 
large values of \au\ do not agree with each other psj . 
Here we shall fill this gap, presenting some results for 
the large- 1 aw I expansion of s Ai m . 

A straightforward generalization of Flammer's method 
can be easily found for general s. Define a new angular 
wavef unction Zi m (u) through p3 | 



(5) 



Sl m (u) = (1 - U 2 ) 2 Zl m (u), 



and change independent variable by defining x = v2cm, 
where c 2 = -(aw) 2 . Substitute this in J2J to get: 



iy2cx — m(m + 1) 



cx 

+ (2c - x 2 )Zi m , xx - 2(rn + s + l)xZi r< 



2msx 



2c + x 
0. 



Zlrn + 

(6) 



When c — > 00, this equation becomes a parabolic cylinder 
function. The arguments presented in |2ll l23l l24j lead to 



s A lm = (2L + l)c + 0(c°) , c -> 00 , 



(7) 



where L is the number of zeros of the angular wavefunc- 
tion inside the domain. One can show that 



L 



l-\m\, |m| > |s|, 
J -|*|, \m\< \s\. 



(8) 



Higher order corrections in the asymptotic expansion can 
be obtained as indicated in 21] . However, we will not 
need them here. We have verified Eq. J7J) numerically, 
solving Leaver's angular continued fraction for s Ai m as 
a function of the complex parameter au. Our numerical 
results (which will be presented in detail elsewhere) are 
in excellent agreement with previous work HHHEa for 
s = 0. For any s they are consistent with equation Q 
when lor <ti luj and |gkj| is large. 



IV. NUMERICAL RESULTS 

To compute the asymptotic QN frequencies of the Kerr 
black hole we use a technique similar to that described 
in [9(. We fix a value of the rotation parameter a. We 
first compute QN frequencies for which \au>\ ~ 1, so that 
formula J7J) is only marginally valid. This procedure is 
consistent with our previous intermediate-damping cal- 
culations: for example, when we include terms up to or- 



der 



-2 : 



in the asymptotic expansion for oAi m provided 



in |21| . our new results for a ~ 0.1 and I — m — 2 match 
the results for the scalar case presented in 0] at overtone 
numbers 20 < N < 30. Then we increase the overtone 
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tency check. We have tried to fit the resulting "universal 
function" , displayed in Fig. |3 by simple polynomials in 
the BH's Hawking temperature T and angular velocity 
fi (and their inverses). None of these fits reproduces our 
numbers with satisfactory accuracy. It is quite likely that 
asymptotic QN frequencies will be given by an implicit 
formula involving the exponential of the Kerr black hole 
temperature, as in the RN case [lUloll3| . 

For any a, the imaginary part cui grows without bound. 
Quite surprisingly, the spacing between modes 8lui is a 
monotonically increasing function of a: it is not simply 
given by 2ttT, as recent calculations and previous nu- 
merical results suggested 0, |2^| . A power fit in a of our 
numerical results yields: 



5lu! = 1/2 + 0.0438a - 0.0356a 2 



(10) 



FIG. 3: Asymptotic real part ujr = 2vo(a) of the I = m = 
2 gravitational and scalar QN frequencies extrapolated from 
numerical data: u)r — * 2w(l/2) ~ 0.21 as a — *■ 1/2. The 
extrapolated frequencies are independent on I and s. 

index N (progressively increasing the inversion index of 
our "decoupled" continued fraction). Finally, we fit our 
numerical results by functional relations of the form: 

ujr(N) =uj r + J£a + ujfa 2 + Jga 3 , (9) 

where a = 1/Vj or a = ^/l/Vj. At variance with the 
non-rotating case @ , fits in powers of 1 /u>i perform bet- 
ter, especially for small and large a. However, both fits 
break down as a — > 0: the values of higher-order fit- 
ting coefficients increase in this limit, so that subdom- 
inant terms become as important as the leading order, 
and the extraction of the asymptotic frequency ujr be- 
comes problematic. The numerical behavior of subdomi- 
nant coefficients supports the expectation (which has not 
yet been verified analytically) that subdominant correc- 
tions are a-dependent. Therefore one has to be careful to 
the order in which the limits N — > oo, a — > are taken 
[Til IT3I fl8| . These observations are consistent with the 
fact that, in the RN case, the zero-charge limit of the 
asymptotic QN frequency spectrum does n ot y ield the 
asymptotic Schwarzschild QN spectrum [ll], LLd] • 

We have extracted asymptotic frequencies using two 
independent numerical codes. For each value of a, we 
found that the extrapolated value of lor is independent 
of s (Fig. [IJ, independent of I and proportional to m 
(Fig. [5J: ujr = mw(a). 

We obtained ujr computing QN frequencies both for 
scalar perturbations (s = 0) and for gravitational per- 
turbations (s = —2). For defmiteness, in both cases we 
picked I — m = 2. The agreement between the extrap- 
olated behaviors of u>r as a function of a is excellent, 
suggesting that both sets of results are typically reliable 
with an error < 1%. Our results are also weakly de- 
pendent on the number of terms used in the asymptotic 
expansion of s A; m : this provides another powerful consis- 



V. CONCLUSIONS 

Motivated by recent suggestions of a link between clas- 
sical black hole oscillations and quantum gravity, we have 
computed for the first time very highly damped QNMs 
of the Kerr black hole. Our calculation was made possi- 
ble by a decoupling of the radial and angular equations, 
carried out using asymptotic expansions of the angular 
separation constant s Ai m for uj r -C u>i and \au)\ 3> 1. 
Our results are very weakly dependent on the number of 
terms used in the asymptotic expansion of s Ai mi and this 
provides a powerful consistency check. We found that: 

(i) The real part of the QN frequencies lor approaches 
a non-zero constant value. This value does not depend 
on the spin s of the perturbing field and on the angular 
index I. It only depends on the rotation parameter a, 
and is proportional to m: 

ujr = mzu(a). (11) 

We determined 137(a) numerically (Fig. I3J), and showed 
that it is not given by simple polynomial functions of 
the black hole temperature T and angular velocity (or 
their inverses). At fixed a, a fit of our numerical data by 
power series in l/u>i and yl/cJ/ suggests that leading- 
order corrections to the asymptotic frequency are prob- 
ably of order 1/cjj. (ii) The imaginary part uij grows 
without bound, the spacing between modes Sivj being a 
monotonically increasing function of a. 

We wish to stress, once again, that the asymptotic fre- 
quency u)r is independent on the spin s of the perturbing 
field: this is consistent with results for highly damped 
QNMs of (charged) RN black holes [HG3. 

By now it is quite clear that the original Hod proposal 
requires some modification. However, the "universality" 
of the asymptotic Kerr behaviour we established in this 
paper is good news. For both charged and rotating black 
holes the asymptotic QNM frequency ujr depends only 
on the black hole geometry, not on the perturbing field. 
If QNMs do indeed play a role in black hole quantization 
this is an essential prerequisite, and it seems to hold. 
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